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Abstract 

Using 22. 7M BB events collected with the BABAR detector at SLAC, we present preliminary 
measurements of the branching fractions for charmless hadronic decays of B mesons into two-body, 
quasi two-body and three-body final states with pions, kaons, and p and oq resonances. In the 
search for exclusive 5° vr+Tr-vrO, we measure ^(5° p=t(770)7r=F) = (28.9 ± 5.4 ± 4.3) x lO'^, 
together with the relative asymmetry ^phys = —0.04 it 0.18 ± 0.02. We also set the upper limits on 
BO /j0(770)7r°, non-resonant 5° ^ vr+vr-vrO, 5° ^ a^(^ T/7r±)7r=F and 5° ^ K^K^. 
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1 Introduction 



Measurement of the rates and CP asymmetries for B decays into the charmless final states can 
be used to constrain the angle a of the unitarity triangle Q. In the case of three body ■k~^'k~tt^ 
decays, such measurements of a would exploit interference between the — > p^tt"^ modes and the 
colour-suppressed -B — > 

In the case of B^ cl^t^^ , the absence of second-class currents, together with the assumption of 
factorization, provide new constraints on CP observables. The kinematics do not allow interference 
between the oppositely-charged resonances in the Dalitz plot as in the 5° /5(770)7r, but in the 
absence of second-class currents might lead to enhanced direct CP violation Q. 

In the case of B^ — > K^K^, the decay rate is expected to be small (10~^ — 10^^) in the Standard 
Model 1^. Final-state rescattering effects can lead to enhancement of the branching fraction and 
the possibility of large strong phases, with correspondingly large CP-violating charge asymmetries 
0, ^. Observation of the K^K^ decay mode would provide important information about the 
strength of final-state rescattering in charmless B decays. 

2 Analysis 

The data sample used consists of 22.74 million BB events, collected at the PEP-II asymmetric 
e^e~ collider at SLAG, with the BABMl detector Q. Hadronic events are selected based on track 
multiplicity and event topology. We use only good quality tracks. Candidate mesons are 
reconstructed from pairs of oppositely-charged tracks that form a well-measured vertex and have 
invariant mass within 3.5 standard deviations (cr) of the nominal mass Q. Candidate photons 
are defined as showers in the electromagnetic calorimeter that have the expected lateral shape and 
are not matched to a track. Candidate vr*^ mesons are reconstructed by combining pairs of photons 
with an invariant mass between 100 and 160MeV/c^; the vr*^ candidates are then kinematically 
fitted with their mass constrained to the nominal vr*^ mass |^]. Pion candidates (except and oq 
daughters) are required to fail kaon selection criteria. 

We reconstruct the decay B aoir in the mode ao rjir, rj — > 77. To be associated with an 
Tj decay a pair of candidate photons is required to satisfy 0.470 < m^-^ < 0.615 GeV/c^ and the rj 
center-of-mass (CM) momentum must be larger than 0.9GeV/c. The pion track and 77 candidate 
form an oq candidate if 0.90 < m^7r < 1.08GeV/c^. 

Candidate B mesons are selected by exploiting the kinematic constraints provided by the T{AS) 
initial state. First we define an energy-substituted mass ruEs, where ^/s/2 is substituted for 
candidate's energy [|. The second variable used is the difference Ai? between the S-candidate 
energy and ^/s/2. For all modes the niES resolution is dominated by the beam energy spread and is 
approximately 2.5MeV/c^, while IS.E resolution is mode dependent and dominated by momentum 
resolution. Candidates are selected in the range 5.2 < mEs{fnEc) < 5.3GeV/c^ and accepted, 
depending on the decay topology, in various AE ranges, restrictive enough to suppress background 
due to other types of B decays. 

The largest source of background is from random combinations of tracks and neutrals produced 
in the e~^e~ — > qq continuum (where q = u,d,s or c). In the CM frame this background typically 
exhibits a two-jet structure. In contrast, the low momentum and pseudo-scalar nature of B mesons 
from T{4S) decays leads to a more spherically symmetric event. This topology difference is exploited 

^In _B° — > ooTT analysis mss is replaced with the energy-constraint mass niEC = •\/s/4 — where ps is obtained 
by applying kinematic constraints to the four-momenta of the B daughters. 



using event-shape quantities. The first variable is the angle 9t between the thrust axes, in the 
CM frame, of the B candidate and the remaining tracks and photons in the event (ROE). We 
require Icos^tI < 0.9. Another quantity used is a Fisher discriminant J^, a linear combination of 
several discriminating variables like the scalar sum of the momenta of the ROE fiowing into nine 
concentric cones centered on the thrust axis of the B candidate, in the CM frame Another 
set of discriminating variables is defined by L^^'""^ = Pi ^ Icos^jp, which are the momentum- 

weighted sums of the cosines of the angles between the ROE charged tracks (L^^'^ ) or neutral clusters 

(n) 

Lj and the thrust axis of the B candidate. In the analysis of — > aQir these variables are used 
in a non-linear (Neural Network) multi-variate analysis. 

Global detection efficiencies, including branching fractions of intermediate states, are listed in 
Table |l[ Appropriate control samples are used to determine efficiencies for vr*^ and Kg reconstruc- 
tion, particle identification, and selection criteria for rriEs and AE. 

Signal yields are determined with either a simple counting analysis, or with a maximum like- 
lihood fit. For the counting analysis, the yield is defined as Ns = Ni — TZN2, where TZ is the 
background fraction of the number of candidates in the signal region to the number in the side- 
band region, A'^i is the number of candidates in the signal region for on-resonance data and A^2 is 
the number of candidates in on-resonance data observed in the side-band region. In the second 
technique, signal yields are determined from an unbinned maximum likelihood fit using tues or 
niEC: A£J, or AW output, 77 mass (where applicable). In each of the fits, the likelihood for 
a given candidate is obtained by summing the product of event yields and probabilities over all 
possible signal and background hypotheses. Monte Carlo simulated data is used to validate the 
assumption that the fit variables are uncorrelated. The parameters of itles^ nT-EC: and J- PDFs 
are determined from data and are cross-checked with Monte Carlo simulation. 
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Figure 1: Separate region of the Dalitz plot are sensitive to different modes. I: B^ — > p^vr^. II: 
^0 ^ p0^o_ jjj. _^ p/±^T. IV: ^ ^/0^o_ V: ^ charged scalar and 7r=F. VI: B^ neutral 
scalar and vr*^. VII: B^ vr+vr^vr'^ at high mass. 

Data for the B'^ — > vr+vr^vr'^ final state can be represented on a Dalitz plot (see Fig. ||). We sub- 
divide the Dalitz plot into distinct regions, each of which chosen to be sensitive to a single resonance 
such as the /?(770), /?(1450) and /o(400 — 1200). The regions are defined using the invariant mass 
of vrvr-pair combinations and the pair helicity angle defined as the angle between the direction of 



one of the pions and the direction of the parent B meson candidate computed in the vrvr-pair rest 
frame. A counting method is used in this analysis. 

There are four decay rates that are of interest for the decay mode P^t^^, defined by 

TpTT = T{B^ /'^TT ) and Tirp = r{B^ P ^r^) together with their CP conjugates fpjr and f ^rp. 
A non-zero value for the asymmetry, given by: 

^ i^piT ~l~ l^Trp) (TpTT ~l" TtTp) 

would signify direct CP violation in at least one of the decays. 



(1) 



3 Results and Systematics 

The results of the fits or the counting method for the various topologies are summarized in Table |l|. 
In those cases where no evidence of signal is found a 90% confidence level upper limit is computed. 
In the case of the counting analysis, we have used the classical method outlined in and we 
have reduced the background estimate and the efficiency by one standard deviation (systematic) 
before making the calculation. In the case of the maximum likelihood analysis, the upper limit on 

the signal yield for mode k is given by the value of for which /g Cduk/ Cmaxduk = 0.90 
where Cmax is the likelihood as a function of n^, maximized with respect to the remaining fit 
parameters. The result is then increased by the total systematic error, and the detection efficiency 
is reduced by its systematic uncertainty in calculating the branching fraction upper limit. The 
statistical significance of a given channel is determined by fixing the yield to zero, repeating the fit 
and recording the change in — 21n£. 

We have made a preliminary measurement of the CP asymmetry in Eq. ^ of ^phys = —0.04 it 
0.18±0.02, which is consistent with zero. Imperfect knowledge of the PDF shapes, of the detection 

^The numerator in Eq. Illis simply the difference of the two direct CP violations (Fp^ — Fp^) and (F^p — f^p) 
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Figure 2: Left plot: tues and AE distributions for B^ p^vr^. Central plot: projection on the 
r/vr invariant mass axis for B^ ao(r/vr)7r analysis. Right plot: ruES distribution for B^ K^K^. 
The curve is the projection of the maximum likelihood fit result. (Projections from likelihood fit 
are obtained after additional requirements on likelihood ratios) 



efficiencies and of the background subtraction (counting method) are the main sources of systematic 
uncertainties on the branching fraction measurements. Uncertainties in the PDF parameterizations 
are estimated either by varying the PDF parameters within la of their measured uncertainties or 
by substituting alternative PDFs from independent control samples and recording the variations 
in the fit results. 

Table 1: Summary of results for detection efficiencies (e), signal yields (Ns), statistical significances 
and measured branching fractions (B). Upper limits are at 90% CL. 



mode 


e{%) 


Ns ± (stat) ± (syst) 


Stat. Sig. (cr) 


^(10-*^) 


p±(770)7rT 


13.5 ± 1.6 


89± 16±6 


5.0 


28.9 ±5.4 ±4.3 


p0(770)7r0 


7.4 ±0.9 


6.1 ±5.8 ±2.8 


1.0 


< 10.6 


7r+7r-7r°(NR) 


7.5 ± 1.0 


-4.2 ±7.3 ±3.8 


N/A 


< 7.3 


ao(r/7r)7r 


32.8 ± 2.4 


i8.iif:J± 1.6 


3.7 


< 11.5 




36.6 ± 4.6 


SAtli ± 3.5 


1.5 


< 7.3 



4 Summary 

We have measured branching fractions for the rare charmless decay /9^(770)7r^ with its 

asymmetry .Aphys and set upper limits on — > p'^(770)7r'^, non-resonant B^ 7r"'"7r~7r'^, B^ — >■ 
«o r?7r±)7rT and 5° ^ K'^K^. 
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